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SUMMARY 


Studies  of  fission-product  gamma  emanation  shortly  after  fission  (1  to 
100  sec)  have  been  initiated.  The  half-lives  that  will  be  established  in 
this  program,  when  used  in  conjunction  with  initial  yields,  will  provide 
the  fission-product  yield  history  for  calculational  models.  The  studies 
undertaken  include  (1)  an  investigation  of  the  transitive  behavior  of  fission 
products  from  approximately  1  sec  after  fission,  (2)  experiments  to 
establish  recoil  range  as  a  parameter  to  use  for  identifying  the  sources  of 
specific  gamma  rays,  and  (3)  design  of  a  Bteady-state  tape  system  for 
catching  and  subsequently  investigating  recoil-fission  products.  This 
latter  study  is  the  result  of  earlier  investigations  and  shows  promise  of 
providing  extensive,  accurate  data  on  fission-product  half-lives. 

Studies  of  diffusion  of  various  nuclideB  in  the- Ca0-Al203-Si02  system,  the 
Na20-Al203-SiC>2  system,  and  Small  Boy  soil  from  the  Nevada  test  site 
have  been  made.  The  use  of  fission  products  from  slow  U-235  fission 
in  situ  has  been  demonstrated  as  a  surveying  method  for  studying  diffusion 
in  various  matrices  doped  with  U-235.  The  results  of  the  studies  are 
consistent  with  the  compensation  law,  particularly  with  the  compensated 
data  developed  previously  in  this  program.  The  liquid-structure  factor 
appears  to  be  important  in  establishing  diffusivities. 

Transpiration  studies  of  the  volatility  of  Cs  from  CaO-Al203-SiC>2  solu¬ 
tions  indicate  that  simple  gas -phase  transport  laws  apply  to  transpired 
Cs.  The  condensed- state  and  vapor-Btate  chemistry  of  this  system  have 
not  been  totally  resolved. 

A  mathematical  description  of  Borption  of  gaseous  components  by  a  parti¬ 
cle,  using  condensation  coefficients  and  condensed- state  diffusion  coeffi¬ 
cients,  is  presented.  This  treatment  suggests  that  condensation  coeffi¬ 
cients  become  important  to  the  process  when  they  are  lower  than  approxi¬ 
mately  2  x  10"k. 

Mass  spectrometric  studies  have  been  made  that  provide  estimates  of  the 
stabilities  of  the  gaseous  species  Sb^^,,  AS4O6,  SbO,  TeO,  and  Te02> 

The  gaseous  species  As02  and  Sb02  were  not  observed,  thus  their 
stabilities  could  not  be  measured. 
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INTRODUCTION 


In  the  past,  fallout  studies  at  Gulf  General  Atomic  have  been  c  oncerned 
with  evaluating  chemical  fractionation  of  fission  products  during  fallout 
formation.  The  studies  presented  in  this  report  are  an  integral  part  of 
this  evaluation,  but  the  area  to  be  studied  has  been  somewhat  changed. 
Henry's  law  studies  are  not,  as  in  former  studies,  the  only  major  con¬ 
cern  of  the  work  presented  here.  In  place  of  Henry's  law  studies,  mure 
effort  has  been  spent  on  the  delineation  of  the  history  of  isotopic  species 
present  in  the  cloud.  This  is  manifest  in  a  program  to  determine  half- 
lives  of  early  fission-product  nuclides.  Although  only  preliminary  studies 
have  been  made  the  results  are  very  promising  and  it  is  anticipated  that 
considerable  knowledge  can  be  obtained  concerning  fission-product  half- 
lives  in  the  time  region  of  1  to  100  sec  after  fission. 

More  concern  is  being  given  the  condensation  process  in  fallout  formation. 
Calculational  studies  of  controlling  factors  for  sorption  of  fission  products 
by  fallout  have  been  made.  A  range  of  importance  of  condensation 
coefficients  is  suggested  herein,  and  experimental  measurements  have 
been  proposed. 

The  diffusion  studies  program  is  considered  to  be  so  important  that  it  has 
been  expanded.  More  work  on  surveying  methods  and  theoretical  studies 
is  presented.  Studies  of  various  matrices  are  important,  and  the  matrix 
problem  is  of  major  concern  in  supplying  data  for  a  diffusion- controlled 
fission-product  sorption  model. 

Thermodynamic  studies  of  various  gaseous  systems  have  been  continued. 
These  mass  spectrometric  studies  are  valuable  in  that  they  provide  the 
basic  data  required  for  estimating  Henry's  law  constants. 
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SHORT-LIVED  FISSION-PRODUCT  STUDIES 


As  a  part  of  this  study  an  investigation  of  fission-product  gamma  spectra 
shortly  after  fission,  as  a  function  of  time  after  fission,  has  been  initiated. 
The  purpose  of  this  study  is  to  delineate  half-lives  of  some  of  the  fission 
products  that  have  not  yet  been  studied  in  detail.  This  information  is 
important  to  a  model  describing  fallout  formation  because  the  chemical 
properties  and  abundances  of  the  isotopes  during  fallout  formation  define 
the  fractionation  of  fission  products  in  fallout.  Given  initial  fission- 
product  yields,  fis sion- product  half-lives  are  used  to  define  the  quantities 
of  each  isotope  present  after  detonation. 

INITIAL  TRANSIENT  STUDIES 

Half-life  studies  may  be  performed  in  several  ways.  A  first  investigation 
using  a  pulsing  TRIGA  reactor  and  a  "rabbit"  tube  has  been  made.  The 
rabbit  tube  is  a  pneumatic  tube  in  which  a  sample  can  be  rapidly  delivered 
to  the  spectrometer-detector  system  from  the  reactor.  The  detector 
system  consists  of  a  4096-channel  analyzer  and  appropriate  computational 
equipment.  The  pulsing-reactor /rabbit-tube  system  had  been  previously 
Bet  up  at  the  TRIGA  facility. 

In  this  experiment,  I  mg  of  natural  U  in  0.2  g  of  pure  C  was  subjected  to 
approximately  2  x  10^  thermal  neutrons /crn^  in  approximately  0.02  sec, 
causing  about  2  x  10^  fissions.  This  sample  was  rabbited  to  a  position 
about  one  foot  from  a  3  cm^  Ge-Li  detector  and  counted  from  the  time  of 
arrival  (approximately  0.6  sec  after  pulsing)  until  2  sec  after  pulsing.  It 
was  recounted  during  succeeding  6-sec  intervals  or  multiples  thereof  for 
about  18  min,  shortening  the  sample  detector  distance  when  appropriate. 
The  resulting  spectra  were  subjected  to  a  computer  search  program  to 
identify  the  various  gamma  rays  according  to  energy  and  to  evaluate  their 
intensity.  Table  1  is  a  listing  of  energies  of  gamma  rays  recovered  by 
the  programs  during  this  study.  To  be  included  in  this  table,  a  peak  must 
have  been  recovered  from  at  least  three  consecutive  spectra. 

None  of  these  peaks  have  been  identified  so  far.  However,  gamma  rays 
from  Cs-138  were  observed  in  the  last  two  spectra.  Some  other  peaks 
have  energies  corresponding  to  reported  values  from  certain  fission 
products,  but  the  effort  afforded  to  this  work  so  far  has  led  to  few  identi¬ 
fications.  This  identification  problem  was  previously  recognized  by 
Gordon,  Haivey,  and  Nakahara. (0 
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TABLE  1 


ENERGIES  OF  RECOVERED  GAMMA  RAYS 


Energy 

(KcV) 

No.  of  Spectra 
Including  This  Peak 

Time  of  First  Recovery 
(sec) 

85 

18 

0.6 

104 

7 

31 

Ill 

7 

258.  5 

127 

5 

381.  5 

132 

10 

0.  6 

146 

10 

2.  5 

160 

1 1 

101 

169 

4 

0.  6 

178 

6 

320 

200 

10 

101 

206 

6 

2.  5 

227 

15 

2.  5 

239 

10 

75.  5 

266 

8 

209 

292 

4 

320 

310 

6 

46 

318 

11 

101 

349 

7 

132.  5 

403 

6 

10 

412 

10 

75.  5 

438 

7 

258.  5 

458 

7 

75.  5 

464 

5 

381.  5 

544 

7 

2.  5 

590 

13 

31 

L  642 

7 

132.  5 

|  658 

3 

664.  5 

721 

8 

164 

789 

4 

443 

830 

5 

31 

912 

6 

258.  5 

968 

9 

2.  5  j 

1195 

6 

207.  5 

1236 

7 

207.5 

Y- 

3 

. . .  "  'n  '  . .  '  i 

1 

/ 

1 

■« 

One  uf  the  major  problems  with  the  -4096- channeL-analy /.ed  data  is  the 
loi-.o  signal-processing  time  (approximately  100  microseconds).  This 
means  that  a  maximum  of  only  10,000  pulses/ sec  can  be  received,  which 

is  a  weak  signal  when  divided  into  the  individual  gamma  peakB.  More 
extensive  data  will  be  required  for  accurate  studies,  presumably  from 
steady-state  system  studies.  Ln  addition  to  providing  more  accuracy  on 
the  presently  recovered  peaks  of  Table  1,  many  more  peaks  will  be 
recovered  using  a  steady-state  system. 

Another  problem  the  experiment  suffers  from  is  considerable  background 
signal.  Bremsstrahlung  is  believed  to  be  a  major  cause  of  this  problem. 

RECOIL- RANGE  STUDIES 


In  studying  fission  products  shortly  after  fission,  we  have  encountered  the 
problem  of  identifying  the  elemental  sources  of  the  gamma  rays.  The 
normal  way  of  attacking  the  species  identification  problem  is  to  perform 
chemistry  to  separate  fission  products  and  to  then  observe  the  decay  of 
the  identified,  separated  elemental  species.  Project  goals  of  obtaining 
defined  spectra  between  1  sec  and  a  few  minutes  after  fission  would  seem 
to  preclude,  at  Least  for  the  shorter  times,  the  chemical  approach.  Thus, 
an  investigation  oi  the  possibility  of  using  the  physical  separation  associ¬ 
ated  with  fission  recoil  has  been  made. 

In  the  first  experiment,  a  thin  fission  source  (dissolved  normal  U)  cast  in 
cellulose  acetate  (approximately  3  microns  thick)  backed  on  a  6  -micron 
Mylar  film  was  used  as  a  source  of  recoiling  fission  products.  This  film 
was  sandwiched  between  other  6-micron  Mylar  films  and  irradiated  in  a  GGA 
TRIGA  reactor.  About  1/2  hr  after  irradiation,  or  as  soon  as  the  films 
could  be  handled  easily,  the  film  package  was  undone  and  the  two  fiLms, 
one,  two,  and  three  layers  ren.  ved  from  the  source  film,  were  counted 
using  the  Ge-Li  detector  at  appropriate  time  intervals.  The  source  and 
film  arrangement  are  shown  in  Fig.  1,  and  the  results  of  these  counts  are 
presented  as  the  ratio  of  counts  in  the  second  and  third  films  corrected 
for  radioactive  decay  losses  (or  gains)  as  a  function  of  mass  number. 

Before  describing  the  results,  it  is  helpful  to  understand  what  might  be 
expected  to  occur  according  to  a  slightly  oversimplified  model  of  this 
experiment.  First,  assuming  the  energy  of  a  fission- recoil  product  is 
deposited  linearly  with  distance  through  an  absorber,  a  particular  fission 
product  originating  at  a  point  source  would  be  stopped  on  the  surface  of  an 
imaginary  sphere  in  the  absorber.  The  radius  of  the  sphere  is  the  recoil 
range  for  the  particular  fission  product.  If  this  sphere  is  intersected  by 
equally  spaced  planes,  the  areas  of  the  portions  of  the  recuil  sphere  be¬ 
tween  these  planes  can  be  compared  (i.  e.  ,  area  is  proportional  to  the 
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recoiled  fission-product  concentration).  The  area  of  the  sphere  between 
two  planes  will  be  the  same  as  between  another  pair  of  equally  spaced 
planes  intersecting  the  sphere.  The  area  between  a  plane  inter  seeling  the 
sphere  and  one  that  does  not  will  be  proportional  to  the  distance  between 
the  last  intersecting  plane  and  a  plane  tangent  to  the  sphere  parallel  to 
the  intersecting  plane.  Of  course,  a  source  plane  is  a  multitude  of  source 
points  with  the  above  properties  so  that  if  the  intersecting  planes  are 
parallel  to  the  source  plane,  the  above  statements  describe  fission-product 
concentrations  for  this  model.  Finite  thickness  of  the  source  plane  alters 
the  picture  predictably.  However,  the  plane  source  description  is  nearly 
true.  In  this  model,  then,  fission-product  recoil  range  may  be  distin¬ 
guished  by  comparing  the  amount  of  a  fission  product  in  the  last  film  that 
catches  the  fission  product  to  the  amount  in  a  film  nearer  the  source 
plane.  Fission-product  recoil-energy  dissipation  is  not  truly  linear  with 
distance;  thus,  precise  range  determinations  are  more  complicated,  but 
the  real  case  parallels  the  above  presentation. 

Fission-product  recoil  range  is  mainly  dependent  on  mass  number  of  a 
fission  product,  its  initial  energy,  and  stopping  properties  of  the  medium 
in  which  recoiling  takes  place.  In  the  nuclear  scission  process  that 
produces  a  pair  of  fission  products,  each  of  the  highly  repelling  nuclei 
attains  a  momentum  of  equal  magnitude  but  opposite  direction.  The 
energy  of  fission  is  divided  mainly  into  kinetic  energy  of  these  two  parti¬ 
cles.  Generally,  the  lighter-fission  products  made  in  fission  are  more 
highly  energetic  and  more  penetrating.  While  this  may  not  be  universally 
correct,  it  is  a  good  first  approximation. 

We  have  observed  recoil-range  variation  with  mass  in  our  experimental 
measurements,  as  shown  in  Fig.  1.  We  believe  that  this  curve  is  . 
descriptive  of  the  particular  fission-product  mass  and  that  with  better 
statistics  we  can  achieve  identification  of  a  fission-product  mass  from  its 
experimental  recoil  range  within  one  or  two  mass  units  for  a  considerable 
range  of  masses.  Gamma  intensity  data  10  hr  after  fission,  presented  in 
this  figure,  fall  on  the  curve  where  the  statistics  are  good.  In  some 
instances  this  relationship  appears  to  fail.  Where  it  does  fail  (as  with 
Xe-135),  we  have  found  such  experimental  problems  as  unresolved  gamma 
rays,  which  arise  from  two  or  more  fission  products  of  significantly 
different  masses.  This  suggests  that  improvements  in  data  processing 
are  required. 

A  more  extensive  study  of  recoil  range  has  been  made.  A  series  of  Ge-LA 
gamma-ray  spectra,  made  42  min  to  30  hr  after  a  60-min  exposure  of  a 
sample  in  a  GGA  TRIGA  reactor,  has  been  produced.  These  spectra  were 
obtained  from  recoil-fission  products  captured  by  stacked  6-micron  Mylar 
films  —  in  particular,  in  the  second-  and  third-removed  Mylar  films  from 
the  U  sample.  In  Fig.  2,  a  second-film-removed  spectrum  taken  between 
26  and  66  min  after  fission,  is  presented.  The  identities  of  most  of 
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Fig.  2.  Spectrum  for  second  film  taken  26  to  66  min  after  fission 
( 1  -hr  irradiation)  (Sheet  1  of  8) 
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Fig.  2.  Spectrum  for  second  film  taken  26  to  66  min  after  fission 
(1-hr  irradiation)  (Sheet  3  of  8) 
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Fi;'.  2.  Spectrum  for  second  film  taken  26  to  66  min  after  fission 
(1-hr  irradiation)  (Sheet  4  of  8) 
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g.  2.  Spectrum  for  second  film  taken  26 
(1-hr  irradiation)  (Sheet  5  of  8) 
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Fig.  2.  Spectrum  for  second  film  taken  26  to  66  min  after  fission 
(1-hr  irradiation)  (Sheet  7  of  8) 
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the  peaks  in  this  spectrum  are  called  out,  although  a  few  of  these  identifi¬ 
cations  are  tenuous.  The  identifications  that  have  been  made  are  in  agree¬ 
ment  with  the  anticipated  recoil  ranges  and  half-lives  of  the  fission 
products  within  the  accuracy  permitted  by  tht  counting  data.  Thu  data 
from  this  experiment  (this  spectrum  and  nine  others)  are  presented  in 
Table  2  as  a  gamma  energy,  a  peak  intensity  ratio  in  the  third  and  second 
films,  an  apparent  half-life,  and  an  isotope  identification  where  an  assign¬ 
ment  has  been  made.  While  there  may  be  a  few  assignment  errors,  there 
has  been  a  great  deal  of  success  using  the  film  intensity  ratio  in  helping  to 
establish  the  assignment.  The  problem  of  proper  unfolding  of  closely 
spaced  gamma  lines  remains, 

As  a  further  demonstration  of  the  recoil-range  mea  urements,  C3/C2 
ratios  for  1-135  at  five  different  times  after  irradiation  were  measured. 
The  results  are  presented  in  Fig.  3.  The  1-135  peak  intensities  were  too 
weak  at  the  first  time  measured  and  are  not  shown.  The  data  at  the 
second  and  fifth  measured  times  arc  marginal  although  shown  on  the 
graph.  The  deviation  bars  are  calculated  standard  deviations,  o  ,  accord¬ 
ing  to  counting  statistics.  Dropping  four  points  because  of  suspected  bias 
(more  should  be  dropped),  60%  of  the  data  fall  within  one  o  of  C3/C2  =  0.338 
and  85%  within  two  tr,  These  values  should  statistically  be  68%  and  95%, 
respectively.  The  very  accurate  data  points  (2.  5%  it  or  lees)  that  are 
necessary  for  the  projected  studies  are  in  good  agreement  with  this  value. 

The  best  data  on  heavy-mass-i'ission  products  are  presented  in  Fig.  4. 
These  data  certainly  suggest  that  recoil  range  can  be  used  for  nuclide 
identifications  in  the  reported  range  provided  the  relationship  between 
the  nuclide  number  and  C3/C2  iH  monotonic  —  a  line  is  used  to  describe 
the  data  in  the  figure.  The  data  define  the  selected  line  within  one  mass 
unit,  suggesting  that  with  good  count  data  a  nuclide's  maos  number  can 
be  identified  from  its  C3/C2  values  in  this  region,  plus  or  minus  one 
atomic  mass  unit. 

STEADY-STATE  SHORT-LIVED  FISSION -PRODUCT  STUDIES 


The  requirements  for  making  an  extensive,  accurate  study  of  short-lived 
fission  products  can  be  enumerated.  First,  one  must  observe  a  sufficient 
number  of  fission -product  decay  events  at  several  times  after  fission  so 
an  accurate  mathematical  description  of  birth  and  decay  of  fission  prod¬ 
ucts  can  be  obtained.  Of  course,  sufficient  energy  resolution  of  these 
gamma  counts  into  a  specific  disintegration  process  must  be  possible. 
Second,  a  method  of  associating  particular  isotopic  disintegrations  with 
each  observed  set  of  count  data  associated  with  a  particular  gamma 
energy  must  be  available.  The  first  requirement  is  best  satisfied  by 
making  steady  state  studies  of  fission-product  decay  using  a  Ge-Li 
detector  and  an  appropriate  multichannel  analyzer  (4096  channels).  Using 
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TABLE  2 


FISSION -PRODUCT  GAMMA-ENERGY/RECOIL-RANGE  CORRELATION 


TABLE  2  (continued) 


Energy 

(KeV) 

c3/c. 

*1/2 

(hr) 

Isotope 

Energy 

(KcVj 

CJC, 

*1/2 

(hr) 

Isotope 

496 

0.  82 

Ru- 103 

266 

0.  9 

10 

Y-93 

488 

0.  1 

2 

253 

Long 

Zr-97 

468 

0.  8 

4.  4 

Ru-105 

249 

0.  26 

7 

Xe-2-135 

41 1 

0.  22 

Short 

Ba- 141 

240 

0.  9 

3 

452 

0.  22 

Short 

Pr- 146 

228 

0.  388 

Long 

Te-132 

447 

0.8 

3 

Y  -92 

220 

0.  3 

6.  8 

1-135 

433 

0.  3 

1.  9 

Te-1-133 

211 

0.  25 

1.3 

Nd- 149 

430 

0.6 

1.  6 

Sr-92 

201 

0.  39 

0.  8 

Te-134 

42  3 

0.  3 

2.  3 

190 

0.  29 

0.  5 

Ba- J  4 1 

417 

0.  38 

7 

1-135 

181 

0.  42 

1 

405 

0.  5 

1.  2 

167 

0.  25 

1.6 

Ba  -  139 

392 

1.  4 

1-134 

156 

2 

385 

1-132 

150 

0.  4 

0.9 

Te-1-131 

364 

0.  5 

Long 

1-131 

142 

0.  85 

Long 

Tc-99 

357 

0.  88 

Short 

129 

0.  5 

Rh- 105 

354 

13 

116 

2 

343 

0.  4 

Short 

105 

0.  1 

318 

0.65 

Rh-2-105 

92 

Nd- 147 

311 

0.5  6 

1.  3 

86 

2 

306 

0.  71 

0.  5 

Tc- 101 

81 

Short 

Xe  - 1 33 

293 

0.21 

Long 

Ce- 143 

76 

2 

277 

0.35 

1.  0 

Te-134 

59 

30 

274 

0.4 

17 

50 

0.  3 

5 

17 


ATOMIC  NUMBER 


Fig.  4.  Fission-product  recoil-range  parameter  C^/C^ 
in  Mylar  vs,  atomic  mass 
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a  steady-state  system  one  wouLd  be  able  to  observe  the  gamma  radiation 
from  fission  products  a  set  delay  time  after  birth  for  a  long  period.  This 
observation  can  be  accomplished  by  allowing  recoil  of  fission  products  to 
separate  fuel  and  fission  products,  catching  the  recoils  in  some  medium 
and  studying  the  decay  of  these  recoiled  fission  products.  The  second 
requirement  can  be  satisfied  by  using  measured  recoil  ranges  associated 
with  the  established  gamma  energies  to  aid  in  identifying  the  fission- 
product  emitting  the  observed  gamma  ray. 

A  system  that  should  make  measuring  steady-state  recoil  range  and  half- 
lives  possible  is  shown  diagrammatically  in  Fig.  5.  Two  thin  Mylar  tapes 
pass  in  front  of  a  recoiling  fission-product  source.  To  determine  the 
recoil  range  associated  with  a  particular  gamma  ray,  the  counting  rates 
in  the  two  films  need  to  be  well  established  at  a  particular  time  after 
fission.  This  will  be  done  by  passing  each  film  containing  recoiled 
fission  products  from  the  steady-state  fission  source  in  front  of  a  detec¬ 
tor  at  a  specific  delay  time  after  fission  and  continuously  counting  more 
or  less  of  the  moving  film  until  statistically  adequate  count  rates  of 
fission  products  can  be  established.  Once  recoil  range  and,  accordingly, 
gamma-ray  identifications  are  made,  one  of  the  recoil  target  films  will 
be  counted  at  various  delay  times  after  fission  to  provide  a  measure  of 
specific  short-lived  isotope  concentration  versus  time.  From  these  data, 
half-lives  of  the  short-lived  fission  product  (1  to  100  sec)  can  be  derived. 
From  these  measured  half-lives  and  initial  yields,  the  quantities  of 
individual  fission  products  as  a  function  of  time  will  be  better  established 
for  fractionation  estimations. 

Equipment  for  these  studies,  involving  either  neutron-beam  tube  fission 
studies  or  Cf-252  spontaneous  fission  studies  or  both,  has  been  designed 
and  will  be  built  and  used  as  part  of  a  new  contract. 

UNFOLDING  GAMMA-RAY  SPECTRA 


The  various  gamma  rays  of  fission  products  shortly  after  fission  are  so 
abundant  that  measured  spectra  lines  exhibit  considerable  overlap,  even 
with  an  energy  resolution  of  2.  8  KeV  (energy  band  width  at  one-half  peak 
intensity).  Most  of  the  data  reported  here  had  resolution  of  that  order, 
as  shown  in  Fig.  2,  The  studies  made  so  far  require  complex  computer 
processing  for  roughly  one-half  of  the  observed  peaks  to  obtain  reliable 
estimates  of  total  counts  under  a  peak.  For  the  purpose  of  obtaining 
reliable  total  counts  from  these  overlapping  spectra,  a  portion  of  our 
effort  has  been  expended  to  establish  a  least-squares  fit  to  a  gamma-ray 
spectrum  in  the  form  of  sums  of  Gaussian  distributions  and  an  exponential 
background  function.  The  system  has  provided  reasonable  description  of 
some  peak  complexes,  but  it  requires  optimization  and  further  testing. 
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The  data-processing  system  presently  used  is  not  entirely  satisfactory 
for  the  work  we  are  trying  to  accomplish,  particularly  since  peak  areas 
accurate  to  better  than  ±  5%  are  essential  to  the  recoil-range  identifica- 

txun  aluuiea, 


DIFFUSION  OF  RADIONUCLIDES  IN  MOLTEN  SILICATES 


The  diffusion  of  fission-product  nuclides  in  molten  silicate  matrices  is  of 
major  importance  in  any  realistic  model  for  fallout  formation.  (2>  3)  We 
have  continued  our  investigations  of  this  phenomenon.  Studies  of  the  dif¬ 
fusion  of  radiosodium  and  radioantimony  in  molten  matrices  from  the 
sodium  disilicate-albite  system  have  been  completed.  An  investigation  of 
the  diffusion  of  radioantimony  in  molten  Nevada  soil  has  been  performed. 
The  simultaneous  diffusion  of  several  radionuclides  in  molten  media  from 
the  CaO-Al203 -Si02  system  has  been  studied  using  two  different  tech¬ 
niques.  The  latter  study  was  made  possible  by  the  availability  of  a  4096- 
channel  gamma  analyzer  and  a  high-resolution  Ge-Li  detector.  The 
experiments  have  helped  to  clarify  the  role  of  silicate  composition  in 
diffusion. 

DIFFUSION  IN  THE  Na^-Al^-SiQ^  SYSTEM 

The  use  of  a  condensed-state  diffusional  model  to  predict  fractionation  in 
fallout  requires  knowledge  of  diffusion  rates  in  fallout  matrix  material. 
Most  of  the  diffusion  studies  performed  at  GGA  have  been  concerned  with 
the  Ca0-Al203-Si02  system.  Methods  of  predicting  diffusivities  have 
been  developed  and  tested  for  this  system.  For  estimating  fractionation 
during  fallout  formation,  it  would  appear  necessary  to  at  least  test  diffu¬ 
sion-coefficient  estimation  methods  on  other  systems.  This  has  been 
partially  accomplished  in  studies  of  the  Na20-Al£03 -SiO^  system. 

For  these  investigations,  four  compositions  from  the  Na^O- AI2O3 -SiC>2 
ternary  system  were  chosen.  These  silicates  were  100%  albite 
(Na20  •  Al^Oj  •  6  SiO^),  100%  sodium  disilicate  (Na^O  •  2  SiO^),  and 
38:62  and  60:40  mixtures  of  these  silicates.  This  particular  system  was 
chosen  for  several  reasons.  Albite,  a  common  mineral,  can  be  consid¬ 
ered  representative  of  a  medium  refractory  soil  for  fallout  simulation, 
and  sodium  disilicate  is  representative  of  a  low  refractory  soil.  Also, 
there  are  data  reported  for  the  leaching  of  radiocesium-diffused  albite.  ^ 
From  the  latter,  and  using  the  results  of  the  present  study  for  Na,  pre¬ 
dictions  of  diffusivities  of  several  other  radionuclides  can  be  made  for  the 
molten  albite  composition.  This  will  be  discussed  in  detail  below. 

Finally,  this  system  was  chosen  to  study  the  problem  of  precrystalline 
structure,  or  configurational  entropy  effect,  in  molten  silicates. 
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The  first  part  of  this  study  was  concerned  with  the  self-diffusion  of  radio- 
sodium  in  Lnti  molten  Iiidli  xcca,  me  »UU1  eta  liiai  wki  e  uaeti  iui  liic  planer 
source  technique(3)  were  prepared  by  activating  the  silicates  in  a  GGA 
TRIGA  reactor  to  yield  Na-24.  Six  experiments  were  completed  in  the 
temperature  range  of  1072'  to  166  1°  K.  These  points  are  superimposed 
on  the  phase  diagram  reported  by  Schairer  and  Bowen(5)  that  is  approxi¬ 
mately  reproduced  in  Fig.  6.  The  first  four  experime.  ts  with  the  38:62 
and  60:40  matrices  were  completed.  These  four  dir<usivity  values  wer» 
then  combined  with  the  results  for  the  diffusion  of  Na  in  sodium  disilic;  te 
that  were  reported  by  Malkin  and  Mogutnov(6)  to  establish  the  compen  a- 
tion  law  for  this  system.(3)  The  latter  was  obtained  by  treating  the  data 
by  the  method  of  least-squares  to  yield 

log1QDo  =  -  4.  33  +  0.  123  E*  ,  (1) 

where  Do  and  E*  are  the  preexponential  factor  (cm^/sec)  and  activation 
energy  (kcal/mole),  respectively,  in  the  Arrhenius  equation  for  diffusion. 
An  experiment  was  then  done  with  sodium  disilicate  to  compare  with  the 
results  of  Malkin  and  Mogutnov.  At  1248°K,  the  oerimental  value  was 
found  to  be  1.  4  x  10-5  cm^/ sec,  which  compareu  well  with  the  value  of 
1.  5  x  10‘5  cm^/sec  calculated  from  the  equation  reported  by  Malkin  and 
Mogutnov.  An  empirical  correlation  was  then  found  between  E*  and  the 
percent  albite  composition,  and  E*  was  predicted  for  the  100%  molten 
albite  composition  by  extrapolation.  Using  Eq.  1,  D0  was  found  for  that 
matrix  and  an  experiment  was  done  at  16bl'K.  The  experimental  value  of 
the  diffusivity  was  found  to  be  1.6  x  10‘5  cmV sec,  which  compared  well 
with  the  predicted  value  of  2.  5  x  10"5  cm^/  sec.  This  result  substantiated 
the  empirical  correlation  mentioned  above.  The  experimental  results  are 
summarized  in  Table  3. 


TABLE  3 


DIFFUSIVITIES  FOR  TRANSPORT  OF  Na-24  IN 
MELTS  FROM  THE  SODIUM  DISILICATE -ALBITE  SYSTEM 


T 

(8  K) 

Composition 
(%  Albite) (a) 

D  x  10b 

(cmV  sec)(b) 

1072 

38 

2.  8 

1248 

38 

6.  9 

1365 

60 

9.8 

1248 

60 

5.  0 

1248 

0 

14. 

1661 

100 

16. 

(a)  The  uncertainty  is  approximately  2% 

(b)  The  uncertainty  is  approximately  50% 
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g.  6.  Phase  diagram  for  the  sodium  disilicate -albite 

system  showing  parameters  of  diffusion  experiments 


The  Arrhenius  parameters  for  this  system  are  summarized  in  Table  4. 


TABLE  4 

ARRHENIUS  PARAMETERS  FOR  DIFFUSION  OF  Na-24  IN 
MELTS  FROM  THE  SODIUM  DISILICATE -ALBITE  SYSTEM 


Composition 
{%  Albite) 

* 

E 

(kcal/  mole) 

lo8IODo 

(D0  in  cm^/ sec). 

0 

11.9(a) 

-  2.  7 -A1 

38 

13.  5  ±  1.  0(b) 

-  2.  67(d) 

60 

19.  5  ±  1.  0(b) 

-  1.93(d) 

100 

32.  3  ±  1.  5(c) 

-  0.  36(C) 

^  From  Malkin  and  Mogutnov^ 

^  Uncertainties  estimated  from  Eq.  1 

^  From  empirical  correlation  and  Eq.  1 

(d)  ^ 

From  Eq.  1 


Using  the  results  obtained  for  the  diffusion  of  Na  in  the  molten  albite 
together  with  the  leaching  data  for  the  radiocesium-diffused  melts  of 
albite  composition  reported  by  Lane,  (4)  it  was  possible  to  predict  diffu- 
sivities  for  several  species  in  this  matrix.  Lane  allowed  Cs  to  diffuse 
into  the  matrix,  at  different  temperatures  above  the  liquidus.  He  then 
leached  the  samples  using  0.  1  N  HC  1 .  From  his  data  E*  for  the  diffusion 
of  Cs  into  the  matrix  is  estimated  to  be  36.  6  kcal/mole.  It  has  been 
shown(3)  that  E*(or  logioDo)  Is  a  linear  function  of  the  reciprocal  ionic 
radius,  r~l, 

*  - 1 

E  -  a  +  br  (2) 

for  the  alkali  metals  and  I  in  those  cases  where  sufficient  data  exist.  For 
the  present  case,  the  coefficients  in  Eq.  2  were  determined  using  the 
values  of  E*  for  Cs  and  Na.  The  values  of  E*  for  I,  Rb,  K,  and  Li  were 
then  calculated.  Equation  1  was  then  used  to  calculate  the  corresponding 
values  of  logloDo*  These  results  are  summarized  in  Table  5. 
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TABLE  5 


EXPERIMENT  A  L  AND  PR.EDICTED  ARRHENIUS  PARAMETERS 
FOR  DIFFUSION  IN  MELTS  OF  ALBITE  COMPOSITION 


Nuclide 

* 

E 

(kcal/molc) 

lo«10°o 

(D  in  cm  /sec) 
o 

I 

37.  8 

0.  32 

Cs 

36.  6 

0.  17 

Rb 

35.  8 

0.  07 

K 

35,  1 

-  0.  01 

Na 

32.  3 

-  0.  36 

Li 

26.  5 

-  1.  07 

It  is  emphasized  that  the  values  reported  in  Table  5  are  estimates  for  all 
of  the  nuclides  except  Na.  The  estimates  are,  however,  based  on  a 
rational  approach  to  the  problem  of  predicting  diffusivities.  It  is  antici- 
pated  that  the  prediction  will  be  tested  experimentally  for  at  least  one 
nuclide. 

A  study  of  the  diffusion  of  radioantimony  in  molten  media  from  the 
sodium-disilicate -albite  system  has  been  completed.  The  temperature 
range  for  the  experiments  was  1123®  to  1687°K.  Experimental  techniques 
were  identical  with  those  described  for  the  diffusion  of  Na  in  this  system. 
The  results  are  listed  in  Table  6. 

The  character  of  the  soil  that  might  be  involved  in  a  near-earth  detonation 
is  critical  to  fractionation  calculations.  If  condensed-state  diffusion  is  a 
limiting  factor  in  the  formation  of  fallout  particles,  and  we  believe  it  is  at 
least  one  of  the  important  factors,  then  orders  of  magnitude  in  the  d  if  - 
fusivity  as  demonstrated  in  the  data  of  Table  6  become  important.  The 
data  are  also  presented  in  Fig.  7  where  the  diffusivity  is  plotted  on  a 
logarithmic  scale  as  a  function  of  the  reciprocal  temperature. 

The  experiment  at  1344'K  using  the  eutectic  matrix  was  performed,  and 
then,  using  the  compensation  law  as  defined  by  the  data  for  the  albite  and 
sodium-disilicate  com, positions ,  the  diffusivity  in  the  eutectic  matrix  at 
1120®K  was  predicted.  The  predicted  value  was  identical  to  the  experi¬ 
mental  value  of  2.  1  x  10-8  cm^/aec  at  1 1  23®  K.  The  result  was  consistent 
with  the  assumption  of  compensation  for  this  system.  The  results  are 
summarized  in  Table  7. 
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I  /T  (  0  K~  X  10?) 


Diffusion  coefficients  for  the  transport  of  Sb-124  in 
liquids  from  the  sodium  disilicate-albite  system  as 
a  function  of  temperature 


TAB  LE  6 


DIFFUSION  COEFFICIENTS  FOR  THE  TRANSPORT  OF  RADIOANTIMONY 
IN  LIQUIDS  FROM  THE  SODIUM  DISILICATE -ALBITE  SYSTEM 


,  Matrix 

T 

("  K)(a) 

D  x  108 

,  2/  ,  (b) 

(cm  /sec) 

Albite 

mm 

55 

Albite 

XU 

25 

Albite 

1859'- 

15 

Albite 

mSSSm-: 

1.  9 

Sodium  Disilicate 

1495 

410 

Sodium  Disilicate 

1407 

170 

Sodium  Disilicate 

1301 

37 

Sodium  Disilicate 

1195 

5.  1 

Eutectic 

13  44 

86 

Eutectic 

1123 

2.  1 

(  21 ) 

The  uncertainty  is  approximately  10°K 
'  '  The  uncertainty  is  approximately  20% 


TABLE  7 

ARRHENIUS  PARAMETERS  FOR  Sb-124  DIFFUSION  IN  MELTS 
FROM  THE  SODIUM  DISILICATE-ALBITE  SYSTEM 


Composition 
(%  Albite) 

* 

E 

(kcal/  mole) 

lo®  10Do 

2  / 

(D  in  cm  / sec) 
o 

0 

5  2.  2 

2.  29 

38 

50.  8 

2.  20 

100 

87.  6 

5.  0b 

In  Fig.  8,  the  energies  of  activation  for  diffusion  of  Na  and  Sb  are  com¬ 
pared  with  the  sodium-disilicate  phase  diagram.  It  is  emphasized  that  all 
of  llie&e  experiments  were  conducted  in  the  liquidus  region.  For  both  Na 
and  Sb  it  appears  that  there  is  at  least  a  qualitative  reflection  of  the  phase 
diagram  by  the  activation  energy.  This  effect,  although  evident,  is  not  as 
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Fig.  8.  Comparison  of  energies  of  activation  for  diffusion  ir. 

molten  media  with  the  sodium  disilicate -albite  phase 
diagram 
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E*  ( KCAL/MOLE ) 


pronounced  as  that  found  for  the  SbOx*CaO-Al2C>3-Si02  system.  (3)  How¬ 
ever,  the  data  show  evidence  of  precrystalline  structure  in  the  liquidus. 
This  effect  has  been  observed  in  several  other  silicate  systems,  for 
example.in  the  Na20-SiC>2  system  by  Dertev  and  Voronkova. ( ' )  Their 
data  for  the  diffusion  of  Na  in  these  molten  silicates  reflect  the  equilib¬ 
rium  phase  diagram;  i.  e.  ,  maxima  for  the  diffusivity  appear  at  two 
eutectics  and  at  the  composition  corresponding  to  NazO  •  2  SiO^.  This 
interesting  problem  deserves  further  attention  and  we  anticipate  conduct¬ 
ing  experiments  to  further  define  it. 

DIFFUSION  IN  MOLTEN  NEVADA  SOIL 


An  investigation  of  the  diffusion  of  radioantimony  in  molten  Nevada  soil 
taken  from  the  Small  Boy  detonation  site  was  completed.  The  purpose  of 
the  study  was  to  make  a  second  comparison  of  Nevada  soil  with  the 
Ca0-Al203-Si02  eutectic,  which  is  the  reference  matrix  used  in  the  GGA 
laboratory. 

Sixty-day  Sb-124  of  high  specific  activity  was  purchased  from  ORNL.  The 
chloride  solution  was  converted  to  a  nitrate  solution  by  several  treatments 
with  nitric  acid  followed  by  evaporation  to  near  dryness.  This  prepara¬ 
tion  was  added  to  powdered  Nevada  glass,  melted,  redivided,  and 
remelted  to  ensure  homogeneity.  The  final  concentration  of  the  redivided 
source  was  roughly  10-2  wt-%  Sb.  As  in  similar  previous  experiments, (3) 
an  effort  was  made  to  obtain  a  minimum  of  four  points  and  an  order-of- 
magnitude  change  in  the  concentration  gradient  for  each  isothermal  exper¬ 
iment,  and  a  minimum  of  four  points  to  define  the  Arrhenius  coefficients 
for  a  particular  system. 

The  Nevada  matrix  used  for  diffusion  measurements  with  radiocesium( 2) 
was  finely  divided  and  heated  at  approximately  1400°C  for  several  days. 
This  was  done  in  an  effort  to  improve  the  vitreous  quality  of  the  matrix, 
since  the  first  preparation  contained  a  few  small  bubbles  and  crystals.  (2) 
The  glassy  character  of  the  silicate  was  improved  by  this  treatment  and, 
following  microscopic  examination,  a  representative  sample  was  analyzed 
using  X-ray  diffraction.  The  diffraction  pattern  was  compared  with  both 
the  patterns  of  the  original  soil  sample  and  the  soil  that  was  preheated  at 
1000°  C.  Although  a  few  lines  were  perceptible  in  the  diffraction  pattern 
of  the  "glass,  "  they  were  faint  and  diffuse.  The  resolution  of  the  camera 
for  this  work  was  roughly  50A.  The  crystallites  responsible  for  the  lines 
could  not  be  ascertained  although  a  search  of  the  ASTM  index  was  made. 

Experiments  with  the  diffusion  of  Sb-124  in  the  Nevada  matrix  were  done 
in  the  temperature  range  of  1218°  to  1766°  K.  The  results  are  listed  in 
Table  8. 
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TABLE  8 


DIFFUSION  COEFFICIENTS  FOR  THE  TRANSPORT  OF 
RADIOANTIMONY  IN  THE  NEVADA  MATRIX 


T 

(°K)(a) 

D  x  108  . 

(cm^/  sec)^^ 

1766 

69 

1695 

32 

1599 

13 

1499 

6.  0 

1416 

5.  4 

1218 

3.  4 

^  The  uncertainty  is  approximately  10“  K 
^  The  uncertainty  is  approximately  20% 


The  data  from  Table  8  are  also  shown  in  Fig.  9  where  the  diffusivity  is 
plotted  on  a  logarithmic  scale  as  a  function  of  the  reciprocal  temperature. 
In  this  figure,  it  is  apparent  that  the  data  fall  into  two  sets.  Microscopic 
examination  of  the  samples  showed  definite  signs  of  devitrification  in  the 
low-temperature  samples  although  the  high-temperature  samples 
remained  glassy.  It  was  concluded  that  the  low-temperature  data  involved 
diffusion  in  increasingly  devitiified  matrices.  The  annealing  periods  for 
the  three  low -temperature  samples  were  16  hr  (1218°K),  3  hr  (1416°K), 
and  1.75  hr  (1499°K).  These  times  are  sufficiently  long  to  allow  devitri¬ 
fication  if  the  temperatures  are  below  the  liquidus  temperature.  If  this 
interpretation  is  correct,  then  the  point  of  intersection  of  the  two 
Arrhenius  curves  in  Fig.  9,  1290°  C,  should  correspond  to  the  approxi¬ 
mate  liquidus  temperature  of  a  "representative"  glassy  Nevada  soil. 

There  is  some  uncertainty  associated  with  the  assumption  of  Arrhenius 
havior  for  the  low-temperature  data,  which  are  summarized  in  the  fol¬ 
lowing  equations: 


Sb/devitrified  Nevada: 

logioD  = 

-  6.  08  -  (1-  70  x  103/T)  . 

(3) 

Sb/ vitreous  Nevada: 

log10D  = 

1.  48  -  (13.  5  x  103/T)  . 

(4) 
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Equation  (4)  is  to  be  compared  with  the  equation  for  transport  of  radio¬ 
antimony  in  the  1 450°  K  eutectic  of  the  CaO -Al-,0, -SiO-,  ternary  system,  O) 

<-  J  c. 


Sb/ 


vitreous  eutectic; 


=  in  -  (14.  1  X  10  /T)  . 


It  appears  that  the  transport  of  Sb  in  the  vitreous  Nevada  matrix  is  simi¬ 
lar  to  that  in  the  eutectic  matrix.  This  result  strengthens  the  statement 
that  the  CaO- AI2O3 -SiO 2  eutectic  matrix  is  a  fair  model  of  vitreous 
Nevada  soil  for  Small  Boy  calculations.^)  Note  that  the  times  involved  in 
cooling  of  an  atomic -device  cloud  are  small  with  respect  to  the  time 
required  to  complete  experiments  in  the  laboratory,  and  devitrification  of 
a  Nevada  glass  would  probably  not  occur.  In  this  sense,  extrapolation  of 
the  Arrhenius  plot  to  lower  temperatures  from  the  high-temperature  data 
is  justified  for  fallout  calculations  if  the  soil  has  been  at  sufficiently  high 
temperature.  It  is  anticipated  that  further  diffusion  studies  will  be  done 
using  the  Nevada  matrix. 

DIFFUSION  IN  THE  CaO-A^O  ^ -SiQ^  SYSTEM 

A  vaporization  technique  has  been  used  to  study  the  simultaneous  diffusion 
of  several  fission-product  nuclides  in  the  1450“  K  eutectic  of  the 
Ca0-Al203*Si02  system  above  the  liquidus.  The  results  of  these  studies 
were  reported  in  a  paper  submitted  for  publication  in  the  Advances  in 
Chemistry  Series.  (®) 

The  diffusion  of  several  fission-product  nuclides  in  other  silicate  glasses 
from  this  system  has  been  studied  in  the  temperature  range  of  1508“  to 
1813°K  using  the  diffusion-couple  method.  Platinum  capillary  crucibles 
were  partly  filled  by  melting  the  matrix  in  place.  The  source  was  then 
melted  in  place  on  top  of  the  matrix.  Sources  were  prepared  by  dissolving 
93%  enriched  UO^  in  the  matrix  at  the  1%  level  followed  by  activation  with 
neutrons  using  the  GGA  TRIGA  reactor.  The  sources  were  allowed  to 
decay  somewhat  after  activation,  prior  to  being  used.  The  capillaries 
were  annealed  at  1508°  ,  1683°  ,  1773°  ,  and  18l3°K  for  known  periods. 
After  annealing,  the  capillaries  were  sectioned  as  in  the  plane  source 
technique.  (3)  The  capillary  length  was  measured  using  micrometer 
calipers  after  each  sectioning,  and  the  abraded  material  together  with  the 
carborundum  paper  was  gamma -analyzed  using  a  Ge-Li  detector  and  a 
4096 -channel  analyzer.  With  this  analyzer,  it  was  possible  to  monitor 
many  individual  fission-product  peaks,  Often  more  than  one  photopeak  of 
a  particular  species  could  be  studied,  thus  providing  a  measure  of  the 
internal  consistency  of  the  data. 
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The  mathematical  model  used  in  these  investigations  is  the  solution  of  the 


J  i - - -a..:  ...at 


r.CCm*  ai>i^u -*iiiuc pcauciu  uinuoivuy 


iur  uiffubion 


from  an  infinite  source  into  an  infinite  sink;  i.  e.  ,  the  lengths  of  the 
source  and  sink  are  large  with  respect  to  the  diffusion  distance.  The 
pertinent  equation  is 


2  C 


C  erfc 
o 


(x  /2  (Dt)1//2  )  , 


(6) 


where  C  is  the  concentration  of  the  diffusing  species,  CQ  is  the  initial 
concentration  in  the  source,  D  is  the  diffusivity  (cm^/sec)  and  t  is  the 
time.  The  axis  of  symmetry  of  the  concentration  is  at  the  original  inter¬ 
face  x  =  0  of  the  source  and  sink  where  2C  =  CD.  The  value  of  D  is 
obtained  by  plotting  values  of  0/Co  against  the  corresponding  values  of  x 
on  probability  graph  paper  and  finding  the  slope  of  the  best  straight  line. 
Experimental  data  conform  well  to  the  model. 

Two  different  matrices  were  used  in  these  studies.  Matrix  A  wt-%  com¬ 
position  was  41  Si02-22  Ai2C>3  -37  CaO,  and  Matrix  B  was  52  Si02-  18 
AI2O3-3O  CaO.  The  experiments  using  Matrix  A  were  done  at  1773'K. 
These  results  are  listed  in  Table  9. 


TABLE  9 

DIFFUSION  COEFFICIENTS  FOR  TRANSPORT  OF 
FISSION -PRODUCT  NUCLIDES  IN  MATRIX  A  AT  1773° K 


Species 

Atomic  Mass 

D  x  107 
(cm  / sec) 

Te 

132 

3.  5 

I 

131 

3.  8 

Tc 

99 

3.  5 

Zr 

95 

5.  1 

Mo 

99 

3.  7 

La 

140 

6.  1 

Ce 

143 

4.  2 

From  the  data  in  Table  9,  it  appears  that  all  of  the  species  diffuse  at 
roughly  the  same  rate  with  an  average  diffusivity  of  4.  2  x  10*?  cm^/sec  at 
1773°  K.  Thus,  the  critical  temperature  T*  is  approximately  177  3"  K  for 
Matrix  A  if  diffusion  for  this  system  is  assumed  to  be  compensated.  (3) 
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Diffusion  in  Matrix  B  was  studied  at  three  different  temperatures,  1313°, 
1683°  ,  and  1508°K-  The  last  temperature  is  below  the  liquidus  tempera¬ 
ture  ol  this  ternary  composition,  it  was  lound  that  by  extrapolating  the 
Arrhenius  line,  obtained  from  the  two  higher  temperatures,  to  1508°K, 
most  of  the  experimental  diffusivitics  at  this  temperature  were  greatly 
underestimated.  Some  devitrification  of  this  matrix  probably  occurred. 
The  data  at  1508°K  are  not  considered  pertinent.  The  data  for  the  two 
higher  temperatures  are  listed  in  Table  10. 


TABLE  10 


DIFFUSION  COEFFICIENTS  FOR  TRANSPORT  OF 
FISSION-PRODUCT  NUCLIDES  IN  MATRIX  B 


Species 

D  x  108 
(cmV  sec) 

1 8  1 3°  K 

1683°  K 

Xe 

9.  5 

2.  7 

Nd 

19 

4.  9 

Ce 

19 

10 

La 

30 

13 

Te 

12 

3.  1 

Ba 

79 

24 

Zr 

7.  5 

2.  5 

Ru 

28 

4.  9 

Although  the  data  in  Table  10  are  preliminary,  they  afford  approximate 
values  of  the  coefficients  in  the  diffusion  equation, 

log10D  =  a  -  (b/T)  .  (7) 

The  values  of  a.  and  b.  are  listed  in  Table  1 1. 
l  l 

The  data  in  Table  11  appear  to  be  compensated,  i.  e.  ,  a  decreases  with 
b  in  a  fairly  regular  fashion.  In  Fig.  10,  the  compensation  correlation  is 
shown  for  these  data.  The  line  drawn  through  the  points  is  a  least- 
squares -fit  to  the  data  with  unit  weighting  assuming  a  linear  correlation. 
This  observation  indicates  that  the  results  of  these  preliminary  experi¬ 
ments  are  internally  consistent  and  that  compensation  prevails  for 
Matrix  B.  It  was  also  calculated  that,  on  a  statistical  basis,  the 
"critical"  temperature  for  this  system  is  1808°K  and  the  "critical"  dif- 
fusivity  is  1.  9  x  lO"?  cm^/sec, 
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13 


17 


Fig.  1 


B  X  10"3  (°K) 


Compensation  correlation  for  the  transport  of 
fission-product  species  in  Matrix  B 
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TABLE  1  1 


COEFFICIENTS  IN  THE  DIFFUSION  EQUATION  FOR 
TRANSPORT  IN  MATRIX  B 


Species 

a 

b  x  10‘3 
(°  K) 

Ru 

3. 29 

17.  8 

Te 

0.  86 

14.  1 

Nd 

0.  78 

13.  6 

Xe 

0.  10 

12.  9 

Ba 

0.  63 

12.  2 

Zr 

-0.  90 

11.  3 

La 

-1.  73 

8.  7 

Ce 

-3.  23 

6.  3 

Although  the  results  reported  here  are  preliminary,  it  is  possible  to 
compare  them  with  previously  reported  results!**)  for  the  1450°K  eutectic 
of  the  CaO -AI2O3 -Si02  system.  In  a  qualitative  sense,  diffusion  in  the 
eutectic  matrix  is  more  rapid  than  in  the  more  refractory  matrices,  A 
and  B.  A  major  purpose  of  this  study  is  to  attempt  a  quantitative  evalua¬ 
tion  of  the  effect  of  composition  upon  the  simultaneous  transport  of 
several  nuclides  in  silicate  systems.  To  accomplish  this  task,  more 
extensive  data  of  this  kind  are  required. 

COMPENSATION  LAW 


The  utility  of  the  compensation  law,  or  compensation  effect,  in  predicting 
the  diffusion  of  fis sion- product  nuclides  in  silicates  has  been  considered 
to  be  important  in  this  laboratory.  (<3*  3)  a  paper  on  this  subject  has  been 
submitted  for  publication.  (9)  The  experimental  data  obtained  during  this 
year  using  the  plane-source  technique!^)  have  been  combined  with  pre¬ 
viously  reported  data(^’  3)  for  diffusion  in  silicates  obtained  in  this 
laboratory.  The  resulting  set  of  28  datum  points  has  been  treated  using 
the  compensation  effect;  the  correlation  between  log[QD0  and  E>,!/4.  575, 
whore  D0  is  the  pr eexponential  coefficient  (cmV sec),  and  E*  is  the 
activation  energy  (kcal/mole)  in  the  Arrhenius  equation.  The  correlation 
is  shown  in  Fig.  1  1.  The  curve  drawn  through  the  points  is  the  least- 
squares-fit 

log  j qDq  -■  ( -4.  37  +  0.  1  10  EV)  ±  0.  60  ,  (8) 

where  the  uncertainty  is  the  standard  error  of  estimate.  The  point  cor¬ 
responding  to  diffusion  of  radioantimony  in  the  devitrified  Nevada  matrix 
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Fig.  11.  Compensation  correlation  for  data  obtained  at  GGA  for 
diffusion  in  silicates.  Solid  circles  refer  to  data 
obtained  this  contract  year 
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was  omitted  from  the  least-squares  treatment  since  it  is  more  than  four 
standard  deviations  from  the  least-squares  Line,  and  transport  can  cer¬ 
tainly  be  expected  to  be  mechanistically  different  from  that  in  glasses. 

The  importance  of  this  correlation  is  that  it  involves  the  diffusion  of  10 
radionuclides  in  glasses  that  range  in  composition  Irom  low -reiractory 
sodium  disiiicate  to  high-refractory  Nevada  soil  and  CaO- A  1  3 -SiO^  The 

ranges  in  the  coefficients  are  roughly  114  kcal/mole  for  the  activation 
energy  and  14  orders  of  magnitude  for  DQ.  In  a  qualitative  sense,  the 
correlation  can  be  used  to  roughly  estimate  diffusion  coefficients  since 
diffusion  in  refractory  matrices  and  large  nuclides  involves  large  values 
of  no  and  E*.  and  diffusion  in  low- refractory  matrices  and  smaLl 
nuclides  involves  small  values.  Quantitatively,  the  correlation  can  be 
used  to  predict  diffusivities  if  an  independent  estimate  of  either  Dc  or  E* 
can  be  made.  This  has  been  discussed  previously.  (3)  This  correlation 
will  be  used  as  the  reference  for  any  further  calculations  done  using  the 
Norman^)  model  and  the  Korts  and  Norman!1^  c.alculational  scheme. 
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HENRY'S  LAW  CONSTANT  MEASUREMENTS 


/  2  1  1  - 1  31 

In  recent  reports,  1  ’  }  Henry's  law  data  have  been  presented  for 

several  radionuclides  in  equilibrium  with  molten  .matrices  from  the 
CaO  -  A 1 2O3 -SiO^  ternary  system.  The  list  of  radionuclides  that  are  amen¬ 
able  to  experimentation  with  this  system,  using  the  transpiration  technique, 
appears  to  be  exhausted.  The  two  major  limitations  are  satisfactory 
volatility  and  lack  of  interaction  with  the  Pt  apparatus.  It  had  been  antic¬ 
ipated  that  a  study  would  be  performed  using  radiotechnetium.  However, 
Tc  was  found  to  react  with,  or  diffuse  into,  Pt  and  these  experiments 
were  discontinued.  Further  Henry's  law  constant  measurements  may  be 
extended  to  other  matrices  using  fission  products  already  found  amenable 
to  study  in  Pt. 

A  study  was  completed  that  evaluated  the  effect  of  gas-phase  diffusion  on 
transpiration  experiments  as  done  in  this  laboratory.  As  a  test  of  the 
transpiration  method  used  in  studying  Cs  Henry's  law  constants  ,  ^  ^  heli¬ 
um  was  used  in  place  of  Ar  carrier-gas.  This  change  was  made  for  both 
the  high-  and  low-concentration  Cs  samples  described  p»  uvious ly .  ^ 

This  study  was  done  at  approximately  1200°C,  Its  purpose  was  to  deter¬ 
mine  if  Cs  was  being  transported  to  the  collector  in  the  gas  phase.  It  is 
possible  that,  at  the  lower  flow  rate,  transport  limits  might  have  been 
achieved  by  diffusion  along  the  walls  of  the  capillary  instead  of  in  the  gas 
phase.  The  reported  lack  of  a  complete  understanding  of  the 
Cs -CaO- A1  ^C^-SiO^  system,^)  enhanced  the  importance  of  this 
exoeriment. 


Gas-phase  diffusivities  can  be  estimated  from  Hir shfelder' s  descrip¬ 
tion,  ( 1 4) 

3/?  i 

0.001858T  '  [(M1  +M2)/M1M2]2 


(9) 


where  Dj  (  2  is  the  interdiffusion  coefficient,  M  is  the  molecular  weight  of 
a  gaseous  species,  ito|  2  an  arithmetic  mean  cross  section,  ilj-j  is  the 
collision  integral  evaluated  at  the  geometric  mean  of  kT/ tx  (k/ (x  is  a 
Leonard- Jones  force  constant  for  x) ,  T  is  the  temperature,  and  P  is  the 
pressure.  A  ratio  Dj  j  t  3  can  be  developed, 
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(10) 


1,  l 


D 


1.3 


J<M!  +  M2>  M3 
»  (M1  + 

,  V  /f  /,V/12  ' 

'  T  ->  I  I  C  1  /  1 


<Ti  +a3 


where  («j  x/k)*/^  is  taken  as  inversely  proportional  to  Qgjj  ,  a  reason¬ 
able  as  sumption  for  values  of  kT/  « ^  larger  than  1 0  .  Above  1 200“ C  this 
assumption  appears  to  be  adequate.  Cross  sections  and  values  of  ex/k 
were  taken  from  Hirshfelder,  and  (r^g,  which  is  not  critical,  was  taken 
as  4.5A.  The  value  of  Dj  t  z/Di  ,  3,  where  1  represents  Cs-1  37,  2  repre¬ 
sents  He,  and  3  represents  Ar,  was  calculated  to  be  2.87.  Experimen¬ 
tally  comparable  Dj  x  ratios  were  obtained;  2.  2  was  obtained  for  the  high 
Cs  concentrations  and  2.7  for  the  low  Cs  concentrations.  These  ratios 
are  in  excellent  agreement  with  the  theoretical  value  within  the  experi¬ 
mental  error.  Thus,  gas-phase  diffusion  is  considered  to  be  the  control¬ 
ling  process  in  the  Cs  transpiration  studies  at  low  gas-flow  rates. 


Those  new  data  do  not  provide  an  explanation  of  the  lack  of  oxygen  or 
water-pressure  sensitivities  described  in  the  last  report.  ^  They  do, 
however,  confirm  the  method  of  analysis  used  to  describe  these  previous¬ 
ly  reported  transpiration  results. 
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THE  ROLE  OF  CONDENSATION 
COEFFICIENTS  IN  FALLOUT  FORMATION 

To  calculate  the  importance  of  condensation  coefficients  in  fallout  forma¬ 
tion,  a  system  has  been  developed  to  evaluate  diffusion  absorption  in 
spherical  particles  where  a  surface  reaction  plays  a  role. 

Let  us  assume  that  the  surface  condensation  rate  is  proportional  to  the 
condensing  species  pressure  P  (atm)  and  reevaporation  is  proportional  to 
its  surface  concentration  Cs  (g/g).  Thus,  for  an  isothermal  system, 


~dt  =  KqcP  -  K'"eCs  1  (H) 

where  C  is  the  average  concentration  of  the  species  in  a  particle,  t  is  the 
time  (sec),  ac  is  the  condensation  coefficient,  ue  is  the  evaporation  coef¬ 
ficient,  and  K  and  K1  are  constants.  To  simplify  the  system,  »c  and  ae 
will  be  taken  as  equal  to  a,  a  constant  independent  of  concentration  or 
pressure.  At  equilibrium 

O  •  §  *  K«P-K'«C<J  .  (12) 


where  Csf  is  the  equilibrium  surface  concentration.  Thus,  since 
P  =  HCsf  where  H  is  the  Henry's  law  constant  (atm/g/g), 

KH  =  K1  ,  (13) 


§ =  KH“  (l-c=) 

For  an  initial  rate  of  absorption  (Cs  -  0), 


and  according  to  the  Langmuir  evaporation  expression 


(14) 


(15) 


44.  3  3 


(4tr  r2) 


P» 


\/MA 


(16) 


43 


where  m  is  the  quantity  condensed  (g),  r  is  the  radius  (cm)  of  the  sphere, 
M  is  the  molecular  weight  (g/mole)  of  the  gaseous  species,  and  T  is  the 
temperature  { ~  K) .  Now, 


/dm\  _  4irr'p  /dC  \ 

\ dt  A "  3  UA  ' 

where  p  is  the  particle  density  (g/cm3)  and 


<«),■ 


132,  99  Pq.  7m/T 
Pr 


17) 


18) 


anil  thus, 


K  = 


_  132.  99  7m/T 


Pr 


(19) 


It  follows,  then,  that  the  rate  of  absorption  of  a  species  by  the  particle  is 


dC  _  132.  99  Ha  i/m/T  /P 


dt 


pr 


«■=.) 


(20) 


(16) 


To  relate  this  information  to  Newman's  __  description  of  diffusion  and 

surface  emission  (Crank  ^ 17^),  the  term  ^  has  to  be  related  to  — , 
Fick's  law  for  a  sphere  is  written  dr 


dC  _ 
dt 


D_  d_  /  2  dC  \ 

2  Or.  (ri  Or.  )  ’ 

ri  1  V  1  / 


[21) 


where  D  is  the  diffusion  coefficient  in  the  condensed  state.  Now 
3  _  r  r 

4tt  r  dC  -  f  a  7  9C,  _  ,  „  f  3  /  2  9C  \ 

— ; —  T r  “  /  4  ir  r .  — dr.  -  4ttD  /  - —  |r.  -  ]  dr 

3  dt  J 0  l  8t  i  JQ  9r.  ^  i  Or.  J  i 


4tt  Dr 


ac 

_ s 

ar 


(22) 
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Therefore, 


ac 

2 

1  Ztr  D  r 

ac: 

g 

3D 

ac 

i 

at 

4  it  r 

8r 

>• 

ar 

Finally, 


(23) 


ac 

S 

dr 


44.  33  Ha  Vm/'T  /p  ^  \  ,  /p  ^  \ 

- Dfi -  U  ■/  "  U  Cs)  * 


(24) 


as  a  description  of  the  surface  condition  where  condensation  and  evapora¬ 
tion  coefficients  are  important.  NewrnanO^)  has  solved  the  diffusion 
problem  for  this  boundary  condition  and  CrankO?)  gives  the  solution  as 


(25) 


where  F  is  the  fractional  absorption  where  the  sphere  starts  with  zero  (or 
uniform)  concentration,  and 

3  cot  6  -  1  -  rh  .  (26) 

r  n  r  n 


The  solutions  to  this  problem  are  presented  by  Crank  for  various  values 
of  rh.  If  rh  is  more  than  10  there  is  little  difference  between  the  solution 
to  Eqs.  25  and  26  and  that  for  a  strictly  diffus ion -controlled  process. 
Condensation  problems,  then,  generally  become  important  for  rh  less 
than  10.  Higher  values  of  rh  may  become  more  important  at  low  absorp¬ 
tions.  However,  as  a  semiquantitative  estimate  of  the  importance  of  con¬ 
densation  coefficients,  if  H  =  1,  T  =  1673,  M  =  104,  p  =  2.  5,  D  ~  1  x  10"^, 
and  r  -  0.  01,  then  the  a  for  which  condensation  problems  begin  (rh  =  10) 
is  2  x  10~^;  that  is,  only  smaller  values  of  a  will  cause  condensation 
problems.  The  values  employed  were  selected  from  the  middle  of  the 
ranges  of  the  quantities  important  for  fallout  calculations.  However, 
these  ranges  are  wide  and,  thus,  condensation  problems  are  undoubtedly 
more  important  than  would  be  suggested  by  ac  ~  Z  x  10*6  limit.  AdamsO^) 
has  used  values  of  a  this  low  and  lower  to  explain  some  experimental  data. 


In  a  more  exact  application  of  Eqs,  25  and  26  to  the  problem  of  fallout, 
difficulty  in  calculation  at  low  absorptions  is  encountered.  For  fallout 
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problems,  important  values  of  F  are  often  below  0.  1.  This  involves 
values  of  n,  r,  D,  t,  and  pn  such  that  the  series  in  Eq.  Z5  is  slowly  con¬ 
verging.  Attempts  were  made  to  transform  this  equation  or  to  simplify 
it  to  render  calculation  simpler.  First  attempts  nave  been  unsucceasiul. 
It  appears  that  a  computer  program  is  required  to  accurately  evaluate 
this  question  throughout  the  entire  range  of  F  values. 

Another  planned  test  of  this  system  is  the  use  of  the  SLIDER  Code,^ 
described  at  the  April  1968  symposium  on  fallout,  to  calculate  the  effects 
of  combinations  of  gas -phase  diffusion,  condensed-state  diffusion,  and 
condensation  rate  on  absorption  by  a  spherical  particle.  Numerical  cal¬ 
culations  of  this  kind  should  be  useful  in  describing  recent  work  by 
Adams!  18), 
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MASS  SPECTROMETRY  OF  OXIDE  SYSTEMS 


The  molecular  character  of  gas eous  - fis  sion  products  and  their  oxides  is 
of  critical  importance  in  describing  fallout  formation.  These  considera¬ 
tions  are  involved  in  many  of  the  mechanisms  of  particle -gas  interaction; 
e.  g.  ,  gas-phase  diffusion,  condensation,  Henry's  law  equilibria, 
condensed-state  diffusion,  and  reevaporation.  During  the  past  year,  the 
method  of  high-temperature  Knudsen-eell  mass  spectrometry  has  been 
used  to  define  gaseous  equilibria  between  metals,  metallic  oxides,  and 
O  for  Sb,  As,  and  Te. 

An  attempt  to  evaluate  the  stabilities  of  MO{g)  and  M02(g)  (M  =  As,  Sb) 
has  been  made.  In  particular,  the  reaction 

M406(g)  rr  4  MO(g)  +  Oz  (27) 

was  studied.  The  metal  oxides  were  dissolved  in  B2O3  to  reduce  their 
activities.  These  melts  were  placed  in  an  Ir  Knudsen  cell  provided  with 
an  oxygen  inlet.  The  cell  was  heated  by  radiation  and  electron  bombard¬ 
ment.  Temperatures  were  measured  by  either  a  Pt-Pt  10%  Rh  thermo¬ 
couple  or  an  optical  pyrometer.  The  molecules  M^jO^g),  MO(g)  and 
were  the  only  species  observed  effusing  from  the  Knudsen  cell  under 
these  conditions.  Appearance  potentials  of  the  oxides  were  found  to  be 
10.  6,  8.  9,  9.  9,  and  9.  0  eV  for  As40^+,  AsO  +  ,  Sb^jO^,  and  SbO+, 
respectively,  with  an  estimated  uncertainty  of  0.  5  eV. 

Isothermal  dependence  of  the  peak  intensity  ratio  M4C>6+:(MO+)^  on  the 
0Z+  signal  showed  that  Reaction  27  was  being  studied. 

Temperature- dependence  measurements  of  the  equilibrium  constants  for 
Reaction  27  were  made.  Four  determinations  were  made  for  each  reac¬ 
tion  with  an  average  of  7  points  each.  Figure  12  and  13  are  typical 
temperature-dependence  plots  for  Reaction  27  for  As  and  Sb. 

Entropies  for  this  reaction  were  calculated  from  equilibrium  constants 
obtained  from  Ag  pressure  calibrations  of  the  systems.  Table  12  is  a 
summary  of  the  heats  and  entropies.  The  data,  including  the  lack  of 
observation  of  MO<>(g),  suggest  that  MO(g)  may  be  the  important  gaseous 
species  for  these  elements  during  fallout  formation.- 


1  r 


Fig.  1  Z .  A  typical  temperature  dependence  of  t.he  equilibrium 
constant  of  Reaction  LI  for  As 
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TAIU.E  I  f 


SUMMARY  OF 

THE  HEATS 

AND  ENTROPIES 

FOR  REACTION 

M 

Temp.  Range 
(°K) 

AH  .  , 

mid  -temp 

( Real/  mole) 

AS  .  , 

mi  d  •  temp 

(eu) 

As 

10^5  -  1175 

5>09  ±  4 

147 

Sb 

U85  -  1415 

3  1  ^  :fc  6 

1  28 

A  study  of  the  reaction 


SbO(g)  .♦  Sb(g)  +  1/iO^ 


(Z8) 


has  also  been  completed. 

A  mixture  of  N.i^OOj  and  Sb^O},  which  had  been  heated  to  its  melting 
point,  was  put  into  an  lr  Knudsen  ceil.  The  cuLl  was  placed  in  the  mass 
spectrometer  and  heated  by  electron  bombardment.  Temperatures  were 
measured  with  an  optical  pyrometer. 

The  ions  0£  +  ,  Sb  +  ,  and  SbO  +  were  obseived  after  considerable  CO  £  was 
evolved.  Appearance  potentials  wer«'  found  to  be  8,  7  and  8.  1  eV  for  Sb  + 
and  SbO  +  ,  respectively.  The  appearance  potential  for  Sb+  is  in  pood 
agreement  with  8.  (>  eV  reported  in  Kiser's  compilation.  The  appear 

aiuc  potential  for  SbO  +  is  t).  9  eV  lower  than  that  found  during  the  study 
of  Reaction  i7.  The  value  of  8.  1  eV  is  believed  to  be  more  reliable 
because  the  SbO  ^  signal  was  a  factor  of  50  higher. 

Four  sets  of  temperature -dependence  measurements  of  the  equilibrium 
constant  for  Reaction  <f8  were  made.  They  yielded  a 

AH°i4i|j  -  58.  1  ±  1.  b  kcal/mole.  The  temperature  range  of  the  measure¬ 
ments  was  1405°  to  15H/°K.  A  silver  calibration  was  done  and  partial 
pressures  were  calculated.  An  entropy  was  calculated  for  Reaction  ^8 
from  the  equilibrium  constant  obtained  from  the  pressure  calibration. 
Table  1  5  is  a  summary  of  the  results. 

A  AH°  [  cq)o  of  4t>7  kcal/mole  is  obtained  for 

Sb  O  (g)  --  4  Sb(g)  +3  0,  (<M) 

4  o  c 


SO 


TABLE  13 


RESULTS  OF  THE  PRESSURE  CALIBRATION  AND 
CALCULATED  ENTROPY  FOR  REACTION  28 


Tern]) 

P 

SbO 

PSb 

PC2 

AS 

(*K) 

(atm) 

(atm) 

(atm) 

(eu) 

1496 

1.  13  x  10" 6 

5.63  x  10 

2.  38  x  id'6 

11.2 

by  combining  the  results  oi  Reactions  27  and  28.  Using  Allj^oo 
60.7  kea!./mol(^)  for  1/2  02  ^  0  and  315  kcal/mole  for  Reaction  27  as 
reported  previously,  values  of  98.8  kcal/mole  for  SbO(g)  and 
831  kcal/mole  for  Sb406(g)  are  obtained  for 

Gaydon^^  reports  a  dissociation  energy,  of  74  ±  9  kcal/mole  at  0UK 
for  SbO(g).  This  value  is  an  extrapolated  spectroscopic  value  in  which 
the  ground  state  of  the  molecule  was  uncertain.  Brewer reports  an 
energy  of  dissociation  to  the  atoms  at  0°K  for  SbijO^g)  of  8H(>  ± 

15  kcal/mole,  a  difference  of  only  4%  from  the  value  found  in  this  study. 

(2  3) 

The  Te-O  system  is  described  elsewhere,. 
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Washington,  D.C.,  20301 

ta  abstract 


Studies  of  fission-product  gamm»  emanation  shortly  after  fission  (1  to  100  sec)  have  been  initiated.  The  half-lives  that  will  be 
established  in  this  program,  when  used  in  conjunction  with  initial  yields,  will  provide  the  fission-product  yield  history  for 
calculational  models.  The  studies  undertaken  include  (1)  an  investigation  of  the  transitive  behavior  of  fission  products  from 
approximately  1  sec  after  fission,  (2)  experiments  to  establish  recoil  range  as  a  parameter  to  use  for  identifying  the  sources  of 
specific  gamma  rays,  and  (3)  design  of  a  steady-state  tape  system  for  catching  and  subsequently  investigating  recoil-fission 
products.  This  latter  study  is  the  result  of  earlier  investigations  and  shows  promise  of  providing  extensive,  accurate  data  on 
fission-product  half-lives. 

Studies  of  diffusion  of  various  nuclides  in  the  CaO-Ai^OySiOj  system,  the  NajO-AljOj-SiOj  system,  and  Small  Boy  soil  from  the 
Nevada  test  site  have  been  made.  The  use  of  fission  products  from  stow  U-235  fission  in  situ  has  been  demonstrated  as  a 
surveying  method  for  studying  diffusion  in  various  matrices  doped  with  U-235.  The  results  of  the  studies  are  consistent  with 
the  compensation  law.  particularly  with  the  compensated  data  developed  previously  in  this  program.  The  liquid-structure  factor 
appears  to  be  important  in  establishing  diffusivities. 

rranspirattor-  studies  of  the  volatility  of  Cs  from  CaO-AljOySiOj  solutions  indicate  that  simple  gavphase  transport  laws  apply 
to  transpired  Cs.  The  condensed-state  and  vapor-state  chemistry  of  this  system  have  not  been  totally  resolved. 

A  mathematical  description  of  sorption  of  gaseous  components  by  a  particle,  using  condensation  coefficients  and  condensed-state 
diffusion  coefficients,  is  presented.  This  treatment  suggests  that  condensation  coefficients  become  important  to  the  process  only 
when  they  are  lower  than  approximately  2  x  \0'b. 


Mass  spectrometric  studies  have  been  made  that  provide  estimates  of  the  stabilities  of  the  gaseous  species  Sb^O^,  As^Og,  SbO,  TeO 
and  Tkij-  The  gaseous  species  AsOj  and  SbOj  weie  not  observed,  thus  their  stabilities  could  not  be  measured. 
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